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I. INTRODUCTION AND OVERVIEW 

Quantum Chromodynamics (QCD) [l], [2] predicts the 
cross sections for the production of objects with large 
transverse momentum (pr) 0] m hadronic collisions. 
Such calculations are performed by factorizing the inter- 
action into a hard scattering process that can be calcu- 
lated perturbatively and a set of soft processes that must 
be described phenomenologically. The high-px jet pro- 
duction cross section is calculated [|| by convolving the 
matrix elements for the scattering of two initial-state par- 
tons (quarks and gluons), with the corresponding parton 
distribution functions (PDF), to produce a partonic final 
state. To predict the momentum spectrum of the final 
particles, additional effects must be considered. The out- 
going partons fragment into jets of hadrons. The beam 
remnants also hadronize and the spectator partons in the 
proton can also interact, leading to multiple parton in- 
teractions (MPI). QCD radiation from the initial- and 
final-state partons occurs, leading to additional jets and 
to an increase in the ambient energy. These effects vary 
with the momentum transfer of the hard parton scat- 
tering. Some of these processes take place at an energy 
scale where the QCD coupling constant is large and per- 
turbation theory cannot be used. They must therefore be 
described using QCD-motivated phcnomenological mod- 
els, implemented in Monte Carlo (MC) event generators. 
In general, these models contain a number of free pa- 
rameters with values that must be obtained by fitting to 
experimental data. 

For a single proton-proton (pp) collision, the Under- 
lying Event (UE) is defined to be any hadronic activity 
not associated with the jets or leptons produced in the 
hard scattering process. In practice, because color fields 
connect all the strongly interacting partons in the proton- 



proton event, no unambiguous assignment of particles to 
the hard scattering partons or UE is possible. Instead, 
distributions that are sensitive to UE modeling are con- 
structed from the tracks that are far from the direction 
of the products of the hard scatter. This direction is ap- 
proximated by the direction of the highest-px (leading) 
object in the event. 

Measurements of observables sensitive to the UE char- 
acteristics in pp collisions were performed at the Teva- 
tron by CDF using jet events at center-of-mass energies 
y/s = 630 GeV and s/s = 1.8 TeV 0] , and using Drell- 
Yan and jet events at y/s = 1.96 TeV @. Underlying 
event observables have been measured in pp collisions 
with y/s = 900 GeV and y/s = 7 TeV from the distribu- 
tion of charged particles in the region transverse to lead- 
ing charged particles and leading charged-particle jets by 
CMS 0,0] and in the regions transverse to and away from 
leading charged particles by ATLAS @ and ALICE 
A complementary analysis by ATLAS studied UE prop- 
erties using both neutral and charged particles fill ]. 

This paper extends previous studies of the UE by mea- 
suring the charged-particle multiplicity and transverse- 
momentum density in the transverse and away re- 
gions with respect to leading charged-particle jets recon- 
structed using the anti-fct [12j | algorithm and varying the 
radius parameter, R, of that algorithm between 0.2 and 
1.0. We can regard the leading charged particle as a 
charged-particle jet with JJrjO. As long as R remains 
below a characteristic value determined by the momen- 
tum transfer (q 2 ) of the hard scatter, the px of charged- 
particle jets formed with larger R are better indicators 
of the hard scatter energy. Such charged-particle jets 
are reconstructed with high efficiency at low jet trans- 
verse momentum, p^', and can therefore be used to study 
the behavior of the UE in the transition region between 
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soft QCD interactions and hard partonic scattering. At 
larger pi^*, charged-particle jets provide a complement to 
calorimeter-based measurements, with results that are 
independent of calorimeter calibrations, selections and 
uncertainties. Variations in the mean values of these UE 
observables with R provide additional information on the 
interplay between the perturbative and non-perturbative 
components of QCD-inspired MC models. 

This paper is organized as follows. The ATLAS de- 
tector is described briefly in Section |TTJ The QCD Monte 
Carlo models used in this analysis are discussed in Sec- 
tion lllll The variables sensitive to UE activity are denned 
in Section IIV1 The event and object selections are pre- 
sented in Section [V] Section PVTl contains a description of 
the analysis. In Section IVlIl the treatment of systematic 
uncertainties is discussed. Fully corrected data distribu- 
tions are presented in Section IVIIII and the results are 
compared to the predictions of several Monte Carlo gen- 
erators. Section IIXI discusses the dependence of the UE 
observables on R. Conclusions are provided in SectionlXl 



II. THE ATLAS DETECTOR 

The ATLAS detector is described in detail in Ref. [l3j . 
The subsystems relevant for this analysis are the inner 
detector (ID) and the trigger system. 

The ID is used to measure the trajectories and mo- 
menta of charged particles. It consists of three detectors: 
a pixel detector, a silicon strip tracker (SCT) and a tran- 
sition radiation straw-tube tracker (TRT). These detec- 
tors are located inside a solenoid that provides a 2 T axial 
magnetic field, and have full coverage in the azimuthal 
angle (j) over the pseudorapidity range \ri\ < 2.5 The ID 
barrel (end-cap) typically has 3 (2 x 3) pixel layers, 4 (2 
x 9) layers of double-sided silicon strip modules, and 73 
(2 x 160) layers of TRT straw tubes. A track traversing 
the barrel typically has 11 silicon hits (3 pixel clusters, 
and 8 strip clusters), and more than 30 straw-tube hits. 
Typical position resolutions are 10, 17 and 130 /im for 
the r-<fi coordinate of the pixel detector, SCT, and TRT 
respectively. The resolution of the second measured co- 
ordinate is 115 /im for the pixel detector and 580 /im for 
the SCT. 

The ATLAS trigger consists of three levels of event se- 
lection: Level-1 (LI), Level-2 (L2), and Event Filter. The 
trigger relevant for this analysis is the single- arm LI min- 
imum bias trigger [li). which uses information from the 
Beam Pickup Timing devices (BPTX) and the Minimum 
Bias Trigger Scintillators (MBTS). The BPTX stations 
are composed of electrostatic button pickup detectors at- 
tached to the beam pipe at ±175 m from the center of 
the ATLAS detector. The coincidence of the BPTX sig- 
nal between the two sides of the detector is used to deter- 
mine when bunches collide in the center of the detector. 
The MBTS consists of 32 scintillation counters located at 
pseudorapidities 2.09 < |?7| < 3.84. The trigger requires 
a single hit above threshold in the MBTS together with 



a BPTX coincidence. 



III. MONTE CARLO SAMPLES 

MC event samples are used to compute detector ac- 
ceptance and reconstruction efficiencies, determine back- 
ground contributions, unfold the measurements for de- 
tector effects, and estimate systematic uncertainties on 
the final results. The samples used here are the same as 
those used in Ref. [Hj]. The baseline Monte Carlo event 
generator used to determine acceptance and efficiencies 
is Pythia 6.4.21 [H| with the ATLAS tune AMBT1, 
which uses the MRST2007LO* PDF 17]. This tune was 
derived using the measured properties of minimum bias 
events jla |. Generated events are simulated using the 
ATLAS detector simulation [l||, which is based on the 
Geant4 toolkit [20]. The simulated events are processed 
using the same software as the data. Several other sim- 
ulated samples are used to assess systematic uncertain- 
ties on the detector response: Pythia 6 using the Peru- 
gia 2010 tune [lj| (CTEQ5L PDF [Hj]), Pythia 6 with 
the ATLAS MC09 tune Q (MRST2007LO* PDF) and 
Pythia 8.145 with tune 4C [H (CTEQ6L1 PDF). 

In the past year, significant work has been done to im- 
prove agreement between the MC generators and LHC 
data by tuning the parameters of the phenomenological 
models used to describe soft QCD processes [24j,|25|. UE 
data, after applying corrections for detector effects, are 
compared to recent tunes of several MC event generators. 
Samples with high statistics were produced using the 
Pythia 6.4.21, Pythia 8.145, and Herwig++ 2.5.1 [26| 
generators. Several different Pythia 6 samples were gen- 
erated with different UE tunes. AUET2B [13 is an AT- 
LAS tune that uses the px-ordered parton shower, inter- 
leaved initial-state radiation (ISR) and MPI which has 
been tuned to UE data from CDF and ATLAS. This 
employs the MRST2007LO** Q3 PDF. The CMS tune 
Zl [25[ is very similar, but is fitted to CMS UE data, and 
uses CTEQ5L PDF. The Perugia 2011 tune Q uses sim- 
ilar settings, with the UE parameters constrained by min- 
imum bias and UE data from CDF and ATLAS. Peru- 
gia 2011 NOCR is a tune of the px-ordered model which 
does not employ color reconnection. The Pythia 8 gen- 
erator includes many new features over Pythia 6, such as 
fully interleaved px-ordered MPI+ISR+FSR evolution, 
and a different model of hard diffraction. The default 
author tune 4C H with CTEQ6L1 PDF [H is used. 
The recently released Herwig++ 2.5.1 version with tune 
UE7-2 used here has angle-ordered parton showers, em- 
ploys MRST2007LO* PDF and has an improved model 
of color reconnection. 



IV. DEFINITION OF VARIABLES 

After reconstructing jets from charged tracks in the 
event, p^' refers to the transverse momentum (px) of the 
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jet with the highest p^. The variation of the UE prop- 
erties with p J r£* is measured in this analysis. In addition, 
the dependence of these properties on R is studied. 

Particles are defined to be in the transverse region if 
their azimuthal angle differs from that of the leading jet 
(\Acj>\ = | ^particle _ ^jet|) by between tt/3 and 2tt/3 radi- 
ans, as shown in Figure [TJ This region is most sensitive to 
the UE. Particles with |A0| > 2ir/3 are defined to be in 
the away region. This region is likely to contain the frag- 
mentation products of the sub-leading parton produced 
in the hard scattering. 

Leading Charged-Particle Jet 
<|> = 




prescaled. The sample represents an integrated luminos- 
ity of 800 /xb _1 after the trigger prescale. The aver- 
age number of collisions per bunch crossing, /i, varied 
throughout the data-taking period, but never exceeded a 
value of n — 0.14. Over half the data were taken with 
/i ^ 0.01. Thus, effects due to the presence of more 
than one collision in the same bunch crossing ( "pile- up" ) 
are minimal. Primary vertex reconstruction |31| is per- 
formed using tracks with pt > 0.4 GeV and |?7| < 2.5. 
A minimum of two tracks is required to form a vertex. 
To further reduce the contributions of pile-up, events are 
rejected if more than one primary vertex is reconstructed. 

Tracks used in the reconstruction of jets and UE ob- 
servables are required to have pr > 0.5 GeV, \n\ < 2.5, 
transverse impact parameter with respect to the primary 
vertex |<io| < 1-5 mm and longitudinal impact parameter 
with respect to the primary vertex |zosin0| < 1-5 mm. 
Only tracks with at least one pixel hit and six SCT hits 
are considered. To minimize the contribution of particles 
produced by secondary interactions in the ID, tracks are 
required to have a hit in the innermost pixel layer if the 
extrapolated track passes through an active portion of 
that layer. 

For each event, jet collections are constructed, corre- 
sponding to the output obtained when the anti-fc t algo- 
rithm is applied to the tracks for five separate values of 
R: 0.2, 0.4, 0.6, 0.8 and 1.0. For each jet collection, the 
leading jet is defined to be the jet with the largest pO^ 1 
satisfying the requirements^* > 4 GeV and 1 77J et | < 1.5. 
This maximum \rf ot \ cut ensures that all tracks associ- 
ated with jets in the fiducial region are within the fully- 
efficient tracking volume. 



FIG. 1. Definition of the transverse and away regions with 
respect to the leading jet. 

Three observables sensitive to UE activity are studied 
in the transverse and away regions: 

• N c h: the number of tracks in the region 

• SpT^ the scalar sum of the transverse momentum 
of the tracks in the region 

• p T : the average pt of the tracks in the region (p T = 

Spx/^ch) 



V. EVENT SELECTION AND 
RECONSTRUCTION 

The events used in this analysis were collected with the 
ATLAS detector at a center-of-mass energy i/s = 7 TeV 
during early 2010. The data sample, event selection and 
reconstruction are identical to those used to measure the 
cross section and fragmentation functions of jets recon- 
structed from tracks (l5j . Events are required to have 
passed the LI minimum bias trigger which was highly 



VI. ANALYSIS PROCEDURE 

The analysis is performed in parallel on each of the 
five jet collections using the following procedure. First, 
the leading jet is selected and events are rejected if there 
is no jet that satisfies the requirements described in the 
previous section. Next, the Acj) of each track with re- 
spect to the leading jet is calculated and the tracks in the 
transverse and away regions are identified. In addition 
to satisfying the selection criteria discussed in Section IVl 
tracks used for the UE measurements are required to pass 
the same pseudorapidity cut as the jets, \rj\ < 1.5. This 
requirement minimizes the contamination in the UE cal- 
culation, from constituent tracks of a leading jet with 
\rf et \ > 1.5. 

The selected tracks are used to calculate the three 
event observables iV c h, S}?t and p T , denoted generically 
as O, in the transverse and away regions. The final re- 
sults presented here are the distributions and mean val- 
ues of these observables for specific ranges (bins) of p^*. 
To allow comparisons with MC generators, the data dis- 
tributions are corrected for detector acceptance, recon- 
struction efficiency, and for bin migration due to track 
and jet momentum-resolution effects. The corrections 
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used in this unfolding procedure are obtained by match- 
ing the jets reconstructed in simulated MC samples with 
those obtained when the anti-fc t algorithm is applied to 
the primary charged particles produced by the genera- 
tor. Primary charged particles are defined as charged 
particles with a mean lifetime t > 0.3 x 10~ 10 s, which 
are produced in the primary collision or from subsequent 
decays of particles with a shorter lifetime. Thus, the 
charged decay products of K s and A particles are not in- 
cluded. A reconstructed jet is considered to be matched 
to a particle-leve l jet if their centers are separated by 
^(Ary) 2 + (A</>) 2 < R/2. 

Jets below the minimum p 3 ^ cut or outside the maxi- 
mum \if Qt \ cut are used in the procedure that corrects for 
the effects of resolution smearing. Thus, looser require- 
ments of p 3 ^ > 1 GeV and \rf ei \ < 2.5 are imposed dur- 
ing the reconstruction. In addition, since corrections for 
migration of jets into the fiducial volume require knowl- 
edge of the population outside that volume, the trans- 
verse momentum of the leading jet outside the fiducial 
volume (p^ xt ) is also determined in each reconstructed 
jet collection for each event and, in MC samples, for each 
collection of particle-level charged-particle jets. 

The events satisfying these requirements are corrected 
back to the primary charged-particle spectra satisfying 
the event-level requirement of at least one anti-/c t jet with 
> 4 GeV and |??> et | < 1.5 reconstructed from charged 
primary particles with px > 0.5 GeV and |?7| < 2.5. Data 
distributions are unfolded using an iterative method [32| 
based on Bayes' theorem, implemented in the RooUn- 
fold [331 ] software package. The procedure requires three 
inputs: a measured input distribution (stored as a multi- 
dimensional histogram), a response matrix (obtained 
from simulated data) that provides a mapping between 
reconstructed objects and those obtained directly from 
the event generator, and an initial choice for the prior 
probability distribution, or prior for short. Each observ- 
able O is stored in a three-dimensional histogram (one 
histogram for each observable and separate histograms 
for the transverse and away regions) where the binning 
variables are p^', O and p™*. To accommodate the de- 
creasing statistics in the data with p^' and the variation 
of the piS and p^ ack resolution with transverse momen- 
tum, these histograms have variable bin width. The re- 
sponse matrix is stored as a six-dimensional histogram 
that specifies the probability that observed values of p^', 
p™* and O are measured for given true values of Pj?, 
p™* and O. This response matrix is not unitary be- 
cause in mapping from generator to reconstruction some 
events and jets are lost due to inefficiencies and some 
are gained due to misreconstruction or migration of true 
objects from outside the fiducial acceptance into the re- 
constructed observables. 

Unfolding the experimental distribution to obtain the 
corrected distribution is done as follows. The response 
matrix, measured data and initial prior are used as in- 
puts to the unfolding algorithm to produce an updated 



distribution, the posterior. This posterior is used as the 
input prior for another iteration of the algorithm, and 
this process is repeated. The inputs to each iteration of 
the unfolding algorithm are the baseline response matrix, 
measured data and the posterior of the previous iteration. 
The optimal number of iterations is determined from MC 
studies. In this analysis, a total of four iterations are per- 
formed for each measured distribution. The initial prior 
is taken to be the prediction of the baseline Monte Carlo 
generator. Systematic uncertainties associated with this 
choice and with the modeling of the response matrix are 
discussed in Section IVIII 

Once the corrected distributions have been obtained 
in bins of p^', the mean value of O for each p^* bin is 
determined from these distributions. Some care must be 
taken to avoid bias when the mean is calculated since 
the output of the unfolding procedure is a histogram and 
the distribution of the population varies across the bin. 
For Ej»t and N c h, the cumulative distribution function of 
the unfolded distribution is calculated and fit to a cubic 
spline and the mean is determined from the results of the 
spline fit. This step reduces the bias between the binned 
and unbinned calculation of the mean from a few percent 
to less than 0.5%. The p T distributions have sufficiently 
fine binning that the bias is below 0.5% without this step. 

VII. SYSTEMATIC UNCERTAINTIES 

A summary of the systematic uncertainties and how 
they affect the measurements is presented in Table HI The 
following sources of systematic uncertainty have been 
considered: 

1 . The track reconstruction efficiency and momentum 
reconstruction uncertainty. 

2. Potential bias arising from the unfolding procedure. 

3. Misidentification of the leading jet due to cases 
where the leading reconstructed jet does not cor- 
respond to the true leading jet. 

4. The uncertainty in the response matrix, which is 
derived using a particular MC sample and thus de- 
pends on the details of the event generator. 

5. Uncertainty in the calculation of the means from 
the distributions of observables due to discretiza- 
tion effects arising from the finite bin-width of the 
unfolded distributions. 

Uncertainties on the tracking efficiency are r\- 
dependent and vary between 2% (for |^ track | < 1.3) and 
7% (for 2.3 < |r/ track | < 2.5 ) 0. The dominant source 
of this uncertainty comes from possible inaccuracies in 
the description of the detector material in the simula- 
tion. The effect of this uncertainty on the measured ob- 
servables is assessed by randomly removing from the data 
a fraction of the tracks consistent with the uncertainty on 
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TABLE I. The systematic uncertainties associated with measurement of the mean values of Ejjt, N c ^ and p T . 



Relative Systematic Uncertainties 




transverse region 


AWAY region 


Source 


S PT (%) 


iVch(%) 


Pt(%) 


Epr(%) 


N ch (%) 


Pt(%) 


Tracking Reconstruction 


2.1-2.5 


2.0-2.3 


0.2-0.3 


2.1-2.6 


2.0-2.3 


0.1-0.2 


Unfolding Procedure 


1.5-6.0 


1.5-4.0 


1.0-4.0 


1.5-5.0 


1.5-4.0 


1.0-4.0 


Leading Jet Misidentification 


< 1.0 


< 1.0 


< 0.5 


< 1.0 


< 1.0 


< 0.5 


Response Matrix 


0.5-2.1 


0.5-1.6 


0.5-1.6 


0.5-2.2 


0.5-1.6 


0.5-1.2 


Discretization Effects 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


Total 


2.7-6.6 


2.6-4.7 


1.3-4.1 


2.7-5.7 


2.6-4.7 


1.3-4.1 



the tracking efficiency and recalculating the observable. 
The resulting uncertainty is then assumed to be sym- 
metric. Uncertainties on the track momentum resolution 
are parameterized as an additional ^-dependent broaden- 
ing of the resolution in track curvature with values that 
vary from 0.4 TeV -1 to 0.9 TeV -1 [Hj]. Systematic un- 
certainties on the tracking performance lead to relative 
uncertainties on the mean values of T,pT and iV c h that 
vary with R from 2.1% (R = 0.2) to 2.6% {R = 1.0). 
Uncertainties on p T are below the percent level for all 
values of R and p 3 ^. 

The performance of the unfolding procedure is stud- 
ied by unfolding the distributions measured in simulated 
MC control samples and comparing them to the known 
generator-level distributions. These closure tests have 
been performed using all the simulated samples described 
in Section [TTTT In each case, before the test is performed, 
the input MC sample is reweighted so that the truth dis- 
tributions of p 3 ^, and of the observable O reproduce the 
unfolded distributions in the data. The resulting devia- 
tions on the mean value of O are lowest (1%-1.5%) at 
low pIS and increase to 4%-6% at large p 3 ^ for all three 
observables. 

At least three effects can result in differences between 
the leading reconstructed jet and the true leading jet. 
One possibility is that a jet with p 3 ^ 1 < 4 GeV or with 
|77 jot | > 1.5 is reconstructed to be inside the acceptance 
("feed- in" jets). Another effect is that due to differences 
in the distributions of true and reconstructed particles 
the anti-A:* algorithm, when applied to the reconstructed 
data, produces jets that do not match any jets obtained 
when the algorithm is applied to true particles ("acci- 
dental jets" ). A third possibility is that a non-leading jet 
is identified as a leading jet due to resolution smearing 
and inefficiencies in the track reconstruction. The un- 
folding procedure intrinsically corrects for migration into 
the fiducial region. Residual uncertainties on the contri- 
bution from misidentification of the leading jet have been 
assessed by reweighting the simulation to reproduce the 
observed distributions of the subleading p^ and the az- 
imuthal angle between the leading and subleading jet, 



and applying the unfolding procedure to the reweighted 
sample. Changes in the resulting output from the un- 
folding differ from the default by < 1%. Studies with 
simulated data indicate that the rate for accidental jets 
is below 0.1% for all values of R and p 3 ^ 1 . Therefore, the 
systematic uncertainty due to accidental jets is judged 
to be negligible. The systematic uncertainty on the frac- 
tion of non-leading jets that are misidentified as the lead- 
ing jet due to uncertainties in the tracking efficiency is 
already included in the tracking systematic uncertainty 
described above. 

Systematic uncertainties due to discretization effects 
(finite bin size) have been studied using simulated data 
by comparing the mean calculated from the binned data 
to the mean calculated from unbinned data. The differ- 
ences are below 0.5% for all observables and all values of 
R. 

The sensitivity to the number of iterations used for 
unfolding is determined by comparing the baseline results 
to those obtained when the number of iterations is varied. 
Changes in the unfolded results are below 0.5% for all the 
observables. 

The response matrix is derived using a particular MC 
sample; therefore, it depends on the details of the event 
generator and tune. The sensitivity of the result to differ- 
ences between the baseline MC sample and the data have 
been studied by comparing the baseline results to those 
obtained when the data are unfolded using an alterna- 
tive response matrix constructed after the MC sample has 
been reweighted to reproduce the raw p 3 ^, and UE observ- 
ables. Relative differences between the baseline measure- 
ments and those obtained with the reweighted response 
matrix are below 0.25% for all three observables. The ef- 
fects of statistical uncertainties in the response matrix are 
evaluated using the bootstrap method [36| to create 50 
statistically independent samples of the reweighted MC 
simulation and repeating the unfolding procedure on the 
data for each sample. The RMS (relative to the baseline) 
of the resulting ensemble of unfolded results is less than 
0.1%. 
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When comparing the mean values of UE observables 
for different values of R, the correlations among the sys- 
tematic uncertainties must be properly treated. Uncer- 
tainties due to track reconstruction efficiency and due to 
discretization effects are fully correlated among the mea- 
surements and thus do not contribute to the systematic 
uncertainty on the ratio of the mean responses measured 
for different R. Uncertainties due to unfolding are par- 
tially correlated. The systematic uncertainties on the 
ratios of UE observables for different values of R are de- 
termined from the deviations from the baseline ratios, 
of the ratios obtained from MC samples where the in- 
put spectra are varied concurrently for all jet collections. 
The uncertainties on the ratios are typically below 1.5% 
except for the highest bin in p^', where the uncertainty 
on the ratio of R = 0.2 to R = 0.6 rises to 6% for Epr in 
the transverse region. 

The same events are used to reconstruct all the jet col- 
lections; therefore, the statistical uncertainties are also 
correlated among the measurements. The statistical un- 
certainties on the ratio of the observables measured with 
one value of R to those measured with a different R are 
obtained by applying a bootstrap method to the data. 



VIII. MEASUREMENTS OF UE 
DISTRIBUTIONS 

The dependence on p!^ of the mean values of the un- 
folded £pt, ^V"ch and p T distributions is shown in Fig- 
ures[2][3]and[4l respectively. The dependence is shown for 
all five values of R for the transverse region. To facilitate 
comparisons with previous measurements, these mean 
values are reported as densities per unit rj-cj). There- 
fore, the mean values of Epx and -^ch are divided by 
2tt = x At], where A</> = 2tt/3 (tt/3 < \Ac/>\ < 2?r/3) 
and Arj = 3 (—1.5 < r\ < 1.5 contributes At] — 3). The 
qualitative behavior of the distributions is the same for 
all five R values. The mean values of Epx rise rapidly 
with p^ for low pi£ and continue to rise slowly for high 

pip*. The mean values of p T exhibit qualitatively similar 
behavior as those of Epx- The mean value of A^h rises 
rapidly with p^ for low p!^ and approaches a plateau 
for high p^ . The systematic differences in the measure- 
ments as R is varied are discussed in Section IIX1 

The unfolded data are compared to several MC gen- 
erators and tunes. In general, the level of agreement 
between the data and MC samples is reasonable, with 
differences below 20% for all observables and all pl^ bins. 
The Pythia 6 Zl sample reproduces the mean values of 
Epx within uncertainites for all p^ bins. This MC sam- 
ple tends to slightly overestimate Ach for p}^ <; 15 GeV, 
and this manifests itself as a slight underestimation of p T 
in the same p^* range. The Pythia 6 AUET2B sample 
tends to slightly underestimate Epx for p^* 20 GeV, 
and overestimate Ep T at higher p^*. Pythia 6 AUET2B 
reproduces N C \ Y for p^ ^ 15 GeV, and overestimates 



V ch for p£* > 15 GeV. This MC sample slightly 
underestimates the mean values of p T in all p^bins. 
Pythia 6 Perugia2011 exhibits reasonable agreement 
with the data EpT distributions ; having tendencies to 
underestimate the data at low p^* and overestimate it 
for p^' <; 15 GeV. This MC sample tends to overesti- 
mate A" c h and slightly underestimate p T . The other tunes 
show somewhat worse agreement for all distributions. 
Pythia 8 tune 4C underestimates all three observables, 
although the agreement is better for p T . Herwig++ tune 
UE7-2 underestimates the Spx and N^. The Pythia 6 
Perugia 2011 NOCR tune is disfavored by the p T data. 

Figures [5] through [7] show the unfolded distri- 
butions (l/Ar ev )dAT ev /dEp T , (l/Ac V )dAe V /dV c h and 
(1 / N ev )dN ev / dp T in the transverse region for three rep- 
resentative values of R and for low (5-6 GeV) and high 
(31-50 GeV) bins of pS ■ Here N cv is the number of 
events in the sample with pJpt > 4 GeV and \if et \ < 1.5. 
The Pythia 6 Zl tune shows differences with respect to 
the data of less than 20%-30% for most distributions and 
p^' bins. The Pythia 6 AUET2B shows discrepancies of 
up to a factor of two for the Epx and N c ^ distributions at 
large values of pp*. For pp* « 5 GeV, the Perugia 2011 
tune undershoots the data for R = 0.2 by 30% — 75% for 
N ch £ 20 and by 20% - 80% for Ep T ^ 20 GeV. The 
differences are smaller for p^* « 30 GeV and for large 
R values Herwig++ shows an excess of events at low EpT 
and low N c ^ for small p^'; at high p^' it does a reason- 
able job of predicting the shape of the distributions but 
understimates the normalization of Spx and N c h distri- 
butions by about 50% . 

The measured distributions in the away region, to- 
gether with MC comparisons, are shown in Appendix lAl 
The EpT, A^h and p T all continue to rise with Pp , as ex- 
pected. This behavior indicates the presence of a second 
jet recoiling against the leading jet. 



IX. DEPENDENCE OF THE UE ON R 

The distributions of the UE observables change with 
the value of R used in the jet reconstruction. The vari- 
ations are summarized in Figure [5J which compares the 
results obtained for the observables in the transverse re- 
gion for R = 0.2, 0.6 and 1.0. For low p 3 ^, the mean 
values of the Spx and A^h densities are largest for the 
smallest value of R, while they are largest for the highest 
value of R at high pip . In contrast, the mean value of p T 
in the transverse region shows little variation with R. 

As noted in Section HI the phenomenological descrip- 
tion of jet production and the UE is complex. This is es- 
pecially true in the low pS region where the distinction 
between the hard scattering process and the soft physics 
associated with the beam remnants is an artifact of the 
model used to parameterize this physics. Nevertheless, 
several general features of jet production are useful for 
interpreting the observed R dependence of the UE ob- 
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(c) R = 0.6, (d) R = 0.8 and (e) R = 1.0. The shaded band shows the combined statistical and systematic uncertainty on the 
data. The data are compared to the predictions obtained with Pythia 6 using the AUET2B and Zl tunes. The bottom insert 
shows the ratio of MC predictions to data for several recent MC tunes. 
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servables. For example, the inclusive jet cross section 
depends on R. Also, leading jets obtained with different 
R parameters will not in general have the same recon- 
structed centroid position, resulting in differences in the 
definitions of the transverse and away regions. Further- 
more, the amount of transverse momentum collected in a 
jet increases with the value of R used for the reconstruc- 
tion. 

Measurements of the inclusive cross section for jets 
reconstructed from tracks show that this cross section 
increases significantly as the radius parameter R is in- 
creased [ljf. Both the Zl and AUET2B tunes of 
Pythia 6 do a remarkably good job of reproducing the 
measurements therein. The number of jets reconstructed 
in the data depends strongly on the value of R used in 
the anti-fci algorithm. For the events passing the trigger 
and primary vertex requirements described in Section fVl 
the fraction of events with at least one reconstructed jet 
with pip > 4 GeV and \rj\ < 1.5 varies from 4.4% to 
19.8% as R is increased from 0.2 to 1.0. These additional 
jets reconstructed for large R largely populate the region 
of low p!Z . The results shown in Figure [5] indicate that 
including these jets reduces the average hadronic activ- 
ity for p J p*^10 GeV, i.e. the additional i?=1.0 jets are 
found in events with less hadronic activity. 

The difference in azimuthal angle between the lead- 
ing jet reconstructed with R — 0.2 and the leading jet 
reconstructed with R = 1.0 has been studied as a func- 
tion of Pp . These studies demonstrate that event re- 
orientation effects due to changes in the reconstructed 
position of the leading jet for different values of R are 
small. For pS ~ 4 GeV, 7% of events are reoriented 
by tt/3 < \A(p\ < 2tt/3. For ^ > 15 GeV the effect is 
much smaller; less than 1% of events are reoriented. 

For cases where jets reconstructed with different R pa- 
rameters are matched, the px of the jet reconstructed 
with larger R will exceed that of the jet reconstructed 
with smaller R, leading to migration of events to bins 
with higher pj?. Many physical processes influence the 
amount of migration as the radius parameter increases. 
These effects include the collection of additional frag- 
mentation particles, the inclusion of additional hadrons 
produced via soft gluon radiation from the final-state 
parton, and the sweeping of particles from the UE into 
the jet cone. Attempts to compensate for the observed 
independence by correcting p^ using the average UE 
energy deposited within the jet cone were unsuccessful. 
This might be due to the fact that there is a correlation 
between the amount of pt migration and the level of UE 
activity in the event because UE activity exhibits long- 
range correlations in rj-ip. Events with higher UE activity 
will exhibit larger px migration as R is increased. Thus, 
the ability of the MC generators to model the variation 
of the UE with R depends not only on how well the gen- 
erator reproduces the mean UE response, but also on 
how well it models the fluctuations in the UE and how 
correlated these fluctuations are spatially. 



X. CONCLUSION 

Observables sensitive to UE activity in events contain- 
ing one or more charged-particle jets produced in pp colli- 
sions at i/i = 7 TeV have been measured. The properties 
of the UE activity have been studied both in the trans- 
verse and away regions. The jets are reconstructed with 
the anti-fct algorithm, with a radius parameter R varying 
between 0.2 and 1.0. Measurements of the evolution of 
the UE activity with R are also presented. Predictions 
from several MC tunes have been compared to the data. 
The predictions from Pythia 6 based on tunes that have 
been determined using LHC data, namely Zl, AUET2B 
and Perugia 2011, show reasonable agreement with the 
data, not only for the mean event activity but also for 
fluctuations in UE activity within events. Other tunes, 
such as Perugia 2011 NOCR, are disfavored by the data. 
The comparison of the predictions from Pythia 8.145 
tune 4C and the Herwig++ 2.5.1 UE7-2 tune to the data 
indicates that additional tuning of UE parameters is nec- 
essary in these cases. The measurements presented here 
provide a testing ground for further development of the 
Monte Carlo models. 
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the predictions of Pythia 6 with the Zl tune; predictions for AUET2B show comparable agreement. 
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Appendix A: UE Distributions in the Away Region 



The dependence on p^* 1 of the mean value of Epr, N ( 



ch 

and p T for all five values of R is shown in Figures [9] 
through 1111 for the away region. The £pr, -^ch and p T 
values all continue to rise with p!^ . The away region 
is expected to have smaller Epx than the towards re- 
gion. This is mainly because the leading jet by definition 
has the highest summed transverse momentum and the 
distribution of the fraction of the momentum carried by 
charged particles is broad. For p 3 ^ > 10 GeV, the mean 
Epx contained in charged particles in the away region is 
typically between 60% (large p 3 ^ and large R) and 100% 
(small p^ and small R) of the leading jet pt- Herwig++ 
overestimates p T for p^^lO GeV. Pythia 6 Perugia 
2011 (without color reconnection) underestimates p T for 
^'£40 GeV. 

Figures [12] through [14] show the unfolded dis- 
tributions (l/N ev )dN ev /dT,p T , (l/N cv )dN cv /dN ch and 
(l/N ev )dN ev /dp T in the away region for three represen- 
tative values of R and for low (5-6 GeV) and high (31- 
50 GeV) bins of p0£ . Here N ev is the number of events in 
the sample with p)^ > 4 GeV and |77 jct | < 1.5. Most of 
the MC models reproduce the shapes of the distributions 
reasonably, with Pythia 6 Zl and Perugia 2011 provid- 
ing the best agreement. Pythia 6 AUET2B predicts 
values of Spx and that are higher than the data at 
large p 3 ^ . 
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FIG. 9. The mean value, per unit 77 — <j>, of SpT in the away region, as a function of (a) for i? = 0.2, (b) for R — 0.4, (c) for 
R — 0.6, (d) for R — 0.8 and (e) for R = 1.0. The shaded band shows the combined statistical and systematic uncertainty on 
the data. The data are compared to the predictions obtained with Pythia 6 using the AUET2B and Zl tunes. The bottom 
insert shows the ratio of MC predictions to data for several recent MC tunes. 
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FIG. 12. Distributions of Ep T in the away region for (a) 7? = 0.2 and 5 GeV < pi£* < 6 GeV, (b) R = 0.2 and 31 GeV < 
p^ < 50 GeV, (c) 7? = 0.6 and 5 GeV < pp* < 6 GeV, (d) i? = 0.6 and 31 GeV < p^ < 50 GeV (e) R = 1.0 and 
5 GeV < pip < 6 GeV and (f) 7? = 1.0 and 31 GeV < Pp* < 50 GeV. The shaded band shows the combined statistical and 
systematic uncertainty on the data. The histograms show the predictions of several MC models. The bottom insert shows the 
ratio of the MC predictions to the data. 
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